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ABSTRACT 


Damped Lyman-a absorbers (DLAs) as a class of quasi-stellar object (QSO) absorption-line systems are currently our most important 
source of detailed information on the cosmic chemical evolution of galaxies. However, the degree to which this information is biased 
by dust remains to be understood. One strategy is to specifically search for QSOs that are reddened by metal-rich and dusty foreground 
absorbers. In this Letter, we present the discovery of a z = 2.60 QSO that is strongly reddened by dust in an intervening extremely 
strong DLA at z = 2.226. This QSO was identified through a novel selection that combined the astrometric measurements from ESA's 
Gaia satellite with extisting optical and near- to mid-infrared photometry. We infer a total neutral atomic-hydrogen column density 
of log N(H r) = 21.95 + 0.15 and a lower limit on the gas-phase metallicity of [Zn/H] > —0.96. This DLA is also remarkable in that 
it exhibits shielded neutral gas that is visible in C1, and it shows tentative detections of CO molecular bands. The spectral energy 
distribution of the QSO is well accounted for by a spectral energy distribution that is normal for a QSO and is reddened by dust from 
a DLA with 10% solar metallicity, a dust extinction of Ay = 0.82 + 0.02 mag, and an extinction curve that resembles that of the Large 
Magellanic Cloud, including the characteristic 2175 À extinction feature. Such QSO absorption-line systems have been shown to be 
very rare in previous surveys, which have mostly revealed sight-lines with low extinction. The present case therefore suggests that 
previous samples have under-represented the fraction of dusty absorbers. Building a complete sample of such systems is required to 


assess the significance of this effect. 
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1. Introduction 


Damped Lyman-a absorbers (DLAs) have been used to study 
the cosmic chemical evolution of galaxies for almost 30 years 
(e.g. Pettini et al. 1990; De Cia et al. 2018). Huge progress has 
been made, and large samples of well-studied systems are now 
available (e.g. Berg et al. 2016, and references therein). From the 
beginning of this line of research, it was clear that despite the 
enormous advantage of the method in obtaining detailed infor- 
mation on the properties of otherwise nearly invisible galaxies, 
the presence of dust would bias the cosmic chemical evolution 
history derived this way. The reason is that DLAs are identi- 
fied in quasi-stellar object (QSO) spectra, and the presence of 
dust inevitably affects the likelihood of identifying QSOs in most 
types of QSO selection (Fall & Pei 1989; Pei et al. 1991, 1999). 

The study of red and/or dust-obscured QSOs also has a long 
and intricate history (e.g. Ostriker & Heisler 1984; Webster et al. 
1995; Benn etal. 1998; Warren et al. 2000; Glikman et al. 
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2006, 2013, 2018; Fynbo etal. 2013; Banerji etal. 2015; 
Hamann et al. 2017). Despite some disagreement among stud- 
ies about the exact fraction of QSOs that is missed due to red- 
dening, it is clear that a substantial number of QSOs are missed 
when they are selected by their optical colours alone. 

A range of studies has attempted to gauge the importance 
of dust bias in DLA samples. Direct measurements of the effect 
of dust reddening have typically found very small amounts of 
excess reddening in samples of QSOs with intervening DLAs 
(Murphy & Bernet 2016, and references therein). These stud- 
ies used to analyse optically selected QSO samples, however. 
Radio-selected samples (which are free from dust bias as long 
as the optical follow-up is 100% complete) have been built 
(Ellison et al. 2001, 2005; Jorgenson et al. 2006), but are small, 
and there is room for a substantial amount of missing metals due 
to dust-bias (Jorgenson et al. 2006; Pontzen & Pettini 2009). 

Our group has specifically aimed at finding exam- 
ples of QSOs that are reddened by intervening dust-rich 
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Table 1. Log of OSIRIS observations. 
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Fig. 1. Contour plots showing the optical and near-infrared point 
sources at b > 60° after the imposed 20 zero proper-motion cut. In 
the upper panel, we display the SDSS g — r vs. UKIDSS J — K colours, 
and in the bottom panel, the WISE W1— W2 vs. UKIDSS J — K colours 
of the sources in the parent sample. Our near- to mid-infrared colour 
criterion is represented by the red dashed line. The optically red targets 
with u—g > 1 and r-z > 0.5 colours are displayed as blue filled circles, 
and GQ 121840832 is marked with a red square. 


DLAs using new QSO-selection techniques (Fynbo et al. 2013; 
Krogager et al. 2015, 2016a). We have found several such sys- 
tems, which shows that dust bias exists and that our meth- 
ods succeed in locating examples of DLAs that have dropped 
out of optically selected QSO samples (Krogager et al. 2016b; 
Fynbo et al. 2017; Heintz et al. 20182). Our ambition is to build 
complete samples of quasars, and within such samples, iden- 
tify examples of very metal-rich sight-lines, which are interest- 
ing for deep spectroscopic follow-up studies. This is challeng- 
ing because of the overlap with cool stars for point sources that 
are red in optical bands. This problem can be alleviated by con- 
structing a complete and (in terms of colour) unbiased sample 
based on astrometric data, specifically from the second Gaia 
data-release (Gaia-DR2; Gaia Collaboration 2018), removing 
sources with significant proper motions that are not extragalactic 
(Heintz et al. 2015, 2018b). 

This Letter is organized as follows. In Sect. 2 we present the 
spectroscopic observations of a newly discovered z = 2.60 QSO 
with celestial coordinates: RA = 12:18:30.10, Dec = 08:32:15.5 
(J2000.0), called GQ 1218+0832, reddened by dust in a fore- 
ground DLA. In Sect. 3 we present our results on the character- 
istics of this DLA, constraining its metallicity and induced QSO 
extinction. We draw our conclusions in Sect. 4. 
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Date (UT) Grism Resolving Exp. time ^ Airmass 
power (s) 

09/12/2018  R1000B 500 2x1000 1.32-1.22 

10/12/2018  R2500V 1680 1x1500  1.92-1.68 

11/12/2018  R2500V 1680 3x1500  1.36-1.12 

17/12/2018 R2500R 1510 2x1350 1.20-1.12 


2. Target selection and observations 


We have built our parent QSO catalogue using a novel com- 
bination (Heintz et al. 2018b) of astrometry from Gaia-DR2 
and photometry from the optical Sloan Digital Sky Survey 
data release 12 (SDSS-DR 12; Eisenstein et al. 2011), the near- 
infrared UKIDSS (Warren et al. 2007), and the mid-infrared all- 
sky WISE mission (Cutri et al. 2013). We confined the sample 
to a region of the sky of Galactic latitude b > 60? to limit 
the contamination by Galactic sources. We only considered tar- 
gets listed as point sources in both SDSS and UKIDSS and then 


only included sources with total proper motions u = 4/2, + u 
consistent with zero within 2c, that is, a signal-to-noise ratio, 
S/N, = H/Her < 2 (following Heintz et al. 2018b, see Fig. 1). 
To be as complete as possible, we initially only imposed the 
zero-proper motion criterion. However, as a significant num- 
ber of stars were still present, we imposed the additional near- 
to mid-infrared criteria shown in the bottom panel of Fig. 1 to 
exclude the large majority of stellar contaminants. Following 
Heintz et al. (20182), we also imposed the optical colour crite- 
riau — g > l and r — z > 0.5 to specifically target dust-reddened 
QSOs at z » 2. 

The object discussed in this Letter, designated 
GQ 1218--0832, was observed as part of a larger spectroscopic 
campaign. It has r = 20.35, g — r = 1.07, and J — K = 0.42 (AB 
mag) and is marked with a red square in Fig. 1. 

Spectroscopic observations were obtained with the OSIRIS 
instrument at the Gran Telescopio Canarias (GTC). We gathered 
spectra using three different grisms in order to better constrain 
the QSO spectral energy distribution, and intervening metal and 
H1 Lye lines. The listed spectral resolution values are based on 
the measured widths of sky-emission lines. The observing log is 
given in Table 1. 

The spectroscopic data were reduced using standard IRAF! 
procedures. The spectra were flux-calibrated using observa- 
tions of spectro-photometric standard stars observed on each 
night. Both the spectra and the photometric data points were 
corrected for Galactic extinction using the dust maps from 
Schlafly & Finkbeiner (2011). All reported wavelengths are in 
vacuum and corrected to the heliocentric rest-frame. To improve 
the absolute flux calibration and correct for slit-loss, we scaled 
the observed spectra to the r-band photometry from SDSS. 


3. Results 


We display the GTC R1000B spectrum in Fig. 2, together with 
the optical and near- to mid-infrared photometry from SDSS, 


! TRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 
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Fig. 2. GTC R1000B spectrum of the z = 2.60 QSO GQ 1218+0832 displayed together with the photometry from SDSS, UKIDSS, and WISE, 
in the u, 2, r, i, z, Y, J, H, Kj, and W1 bands (Warren et al. 2007; Wright et al. 2010; Eisenstein et al. 2011). Overplotted is the composite QSO 
spectrum from Selsing et al. (2016; in blue) and the same composite reddened by Ay = 0.82 + 0.02 mag assuming LMC-type dust at the redshift 
of the DLA (in red). Models with and without the 2175 À dust extinction feature are shown. The inset shows a zoom-in on the blue part of the 
spectrum around the H 1 Lya line from the z = 2.2261 absorber. The red dashed line is the modelled QSO continuum. Also overplotted is the H 1 
DLA Voigt-profile fit, which is included in the dust-reddened QSO composite spectrum, with N(H 1) = 102! 9950.5 cm-?. 


UKIDSS, and WISE. We identify an extremely strong DLA 
(hereafter ES-DLA; Noterdaeme et al. 2014) at a redshift of 
z = 2.2261 based on a range of associated metal lines (see 
Table 2). 


3.1. Extinction 


The spectrum of the background QSO is inconspicuous except 
for its reddening. Because the QSO looks "normal", it is well 
suited for determining the extinction curve of dust in the inter- 
vening ES-DLA. 

To derive the amount of extinction, we used the composite 
QSO template from Selsing et al. (2016), including a range of 
different extinction models. We followed the parametrization of 
Gordon et al. (2003) and fit the extinguished spectrum assuming 
a dust composition similar to the Small and Large Magellanic 
Clouds (SMC and LMC). We find the best fit to be an LMC-type 
extinction curve with a visual extinction of Ay = 0.82+0.02 mag, 
including the 2175 À dust extinction feature. We measure a 
bump strength of Apump = 0.38 mag (following the definition in 
Gordon et al. 2003). 


3.2. Absorption lines 


In the inset of Fig. 2 we show a zoom-in on the damped Lyman-a@ 
absorption line. To derive the H1 column density of the absorber, 
we included a DLA in the dust-reddened QSO model (using the 
approximation by Tepper-Garcia 2006). From the best-fit model, 
we find log N(H 1) = 21.95 + 0.15. 

We also detect a range of metal lines associated with the 
DLA, as well as another intervening system at z = 2.247 seen 
in C iv and Sirv. In Table 2 we list the equivalent width (EW) 
and redshift measurements drawn from the R2500V and R2500R 
spectra. The average redshift of the system inferred from the 
low-ionization lines listed in the table is z = 2.2261 + 0.0004. 


We fitted Voigt profiles to the low-ionisation metal 
lines and derived the following lower limits on the abun- 
dances: [Zn/H] > —0.96, [S1/H] » -1.23, [Fe/H] » —1.89, and 
[Cr/H] > —1.64, using the solar photospheric abundances from 
Asplund et al. (2009). Because of the limited resolution of the 
spectra, these abundances should be considered cautiously and 
only treated as lower limits due to saturation. Taken at face value, 
our measurements indicate an overall metallicity of at least 1096 
solar and significant depletion of refractory elements, for exam- 
ple, [Fe/Zn] ^ -1. 

We detect clear absorption signatures of C1, but because the 
spectral resolution is limited, the transitions from the ground 
state and fine-structure energy levels are blended. Remarkably, 
there are consistent features also at the expected locations of CO 
molecular bands (see Fig. 3). A deeper spectrum with higher res- 
olution is required to confirm the presence of these lines. 


4. Discussion and conclusions 


In Fig. 4 we compare the metallicity and Hr column den- 
sity of the system to the general population of QSO-DLAs 
from De Cia et al. (2018), gamma-ray burst (GRB)-DLAs from 
Bolmer et al. (2019), the extreme case of GRB 080607 from 
Prochaska et al. (2009), and the compiled sample of dusty (Ay > 
0.1 mag) DLAs from Heintz et al. (2018a), including the recent 
discovery by Ranjan et al. (2018). Compared to other QSO- 
DLAs, the absorption system towards GQ 1218+0832 is remark- 
able in several ways. Firstly, such a high H1 column density is 
very rare in the general population of QSO absorbers (Bird et al. 
2017). Even within the population of extremely strong QSO- 
DLAs, this system has higher-than-average metal-line strengths, 
depletion factors, and H1 column density (Noterdaeme et al. 
2014). Its H1 column density is more similar to the column 
densities found in GRB-DLAs (see also Tanvir et al. 2019). 
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Table 2. Absorption lines from the z = 2.2261 + 0.0004 DLA detected 
in the GTC R2500V and R2500R spectra. 


Transition EW, Redshift 
(À) 
Sirv 2 1393 0.41 +0.11 2.2262 
Siv2 1402 0.53+0.09 22262 
CO 41447 0.20+0.11 2.2265 
CO 141477 0.09 +0.06 2.2265 
CO 41509 0.09+0.07 2.2265 
Sin41526 1.2140.07 22257 
C1v 41548 0.48 +0.05 2.2266 
C1v 41550 0.33 +0.05 2.2255 
C1v 41548 0.46+0.04 2.2474 
C12 1560 (9 0.54+0.04 2.2265 
Fe 1n 1608 1.10 +0.06 2.2257 
Ferd 1611 0.1140.04 2.2256 
C1A 1656 0.58+0.06 2.2265 
Aln 21670 1.43+0.05 2.2260 
Nim A 1709 0.15+0.03 2.2259 
Niwd1741 0.22+0.03 2.2260 
NimnA1751 0.19+0.04 2.2263 
Sim A 1808 0.67+0.04 2.2260 
Alma 1854 0.35+0.05 2.2262 
Alma 1862 0.20+0.05 2.2266 
Tin A1910 0.13+0.04 2.2261 
Zn/CruA2026 0.51+0.04 2.2268 
Cr 122056 0.23+0.04 2.2259 
Zn/CruA2062 0.45+0.04 2.2257 
Cr 1122066 0.142 0.04 2.2265 
Fer 42249 0.170.004 2.2261 
Fen 42260 0.36+0.04 2.2265 
Fem 42344 1.93+0.05 2.2261 


Notes. Blended with C iv 4 1550 at Zaps = 2.2474. 
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Fig. 3. GTC R2500V spectrum showing tentative detection of CO A-X 
bands. In the bottom panel, the geometric mean of the profiles is shown, 
revealing the line-profile asymmetry. 
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This DLA is also unusual through the presence of C1 and 
possibly CO absorption. The incidence of C1 absorbers with 
W,(A1560) > 0.4A is about 100 times less frequent than that 
of QSO-DLAs (Ledoux et al. 2015). The presence of strong C1 
lines is an indicator of high metallicity and diffuse molecular 
gas, and Noterdaeme et al. (2018) found that all such systems 
present H and a significant fraction of them exhibit CO at a 
detectable level. This suggests that CO lines are likely to be 
detected in this system as well, in agreement with the presence 
of consistent dips in the spectrum. This system provides further 
evidence that ES-DLAs (with N(H1)> 10?!7 cm"?) lie in the 
transition regime where H1 is being converted into H5 (Schaye 
2001; Noterdaeme et al. 2015; Balashev & Noterdaeme 2018). 

This case, in terms of extinction and molecular content, is the 
most remarkable example so far of a QSO that was missed by 
traditional optical-colour selections due to dust in an interven- 
ing galaxy, as predicted by Fall & Pei (1989). How many quasars 
would be missed if they had a DLA in their line of sight with the 
same level of extinction as in the present case? The answer to this 
question is complex, as the effect of dust depends on the redshift 
of the quasar, the redshift of the DLA, and the colour distribution 
of stellar populations. In this survey, we find that the fraction of 
candidate quasars with the same red optical colours as the system 
presented here is ~1%, compared to the full set of quasars that is 
identified based on the astrometric and near- to mid-infrared selec- 
tion criteria detailed in Sect. 2. Krogager et al. (2019) have carried 
out a comprehensive analysis of the effect of dust in DLAs and 
found that while the fraction of quasars dropping out of selections 
is small, the effect on the derived mass density of neutral hydro- 
gen and metals is much larger, up to a factor of 2 and 5 at z ~ 2.2, 
respectively (see also Pontzen & Pettini 2009). 
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There is little doubt that more such systems are yet to be 
found. An important objective of future work will be to identify 
them. Another major endeavour will be to use the information 
of their incidence rate to build up a more complete picture of 
cosmic chemical evolution. 
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